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ABSTRACT 

The existence of a dark energy component has usually been invoked as the most plausible way 
to explain the recent observational results. However, it is also well known that effects arising from 
new physics (e.g., extra dimensions) can mimic the gravitational effects of a dark energy through a 
modification of the Friedmann equation. In this paper we investigate some observational consequences 
of a flat, matter dominated and accelerating/decelerating scenario in which this modification is given 
by — g{pm,n,q) where g{pm,n,q) is a new function of the energy density pm, the so-called 
generalized Cardassian models. We mainly focus our attention on the constraints from the recent 
Cosmic All Sky Survey (CLASS) lensing sample on the parameters n and q that fully characterize 
the models. We show that, for a large interval of the q — n parametric space, these models are 
in agreement with the current gravitational lenses data. The influence of these parameters on the 
acceleration redshift, i.e., the redshift at which the universe begins to accelerate, and on the age of 
the universe at high-redshift is also discussed. 

Subject headings: Cosmology: theory — dark matter — distance scale 
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1. INTRODUCTION 

Over the last years, a considerable number of high 
quality observational data have transformed radically 
the field of cosmology. Results from distance measure- 
ments of type la supernovae (SNc la) (Perlmutter et al. 
1999; Riess et al. 1998; 2004) combined with Cosmic 
Microwave Background (CMB) observations from bal- 
lon (de Bernardis et al. 2000; O'Dwyer et al. 2003), 
ground (Mason et al. 2003) and satelite experiments 
(Spergel et al. 2003) and with dynamical estimates of the 
quantity of matter in the universe (Calberg et al. 1996; 
Dekel et al. 1997) seem to indicate that the simple pic- 
ture provided by the standard cold dark matter scenario 
(SCDM) is not enough. These observations are usually 
explained by introducing a new hypothetical energy com- 
ponent with negative pressure, the so-called dark en- 
ergy or quintessence (For recent reviews on this topic, 
see Sahni & Starobinsky, 2000; Peebles & Ratra 2003; 
Padmanabhan, 2003; Lima, 2004 and references therein). 
Besides its consequences on fundamental physics, if con- 
firmed, the existence of this dark component would also 
provide a definitive piece of information connecting the 
inflationary flatness prediction with astronomical data. 

On the other hand, it is also well known that not less 
exotic mechanisms like, e.g., geometrical effects from ex- 
tra dimensions may be capable of explaning such observa- 
tional results. The basic idea behind these "braneworld 
cosmologies" is that our 4-dimensional Universe would be 
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a surface or a brane embedded into a higher dimensional 
bulk space-time to which gravity could spread (Hofava 
& Witten, 1996a; 1996b; RandaU L. & Sundrum, 1999; 
Dvali, Gabadadze & Porrati, 2000; Perivolaropoulos & 
Sourdis 2002; Maia et al., 2004). In some of these sce- 
narios the observed acceleration of the Universe can be 
explained (without dark energy) from the fact that the 
bulk gravity sees its own curvature term on the brane 
acting as a negative-pressure dark component which ac- 
celerates the Universe. A natural conclusion from these 
and other similar studies is that dark energy, or rather, 
the gravitational effects of a dark energy could actually 
be achieved from a modification of the Friedmann equa- 
tion arising from new physics. 

Following this reasoning, several authors have re- 
cently investigated cosmologies with a modified Fried- 
mann equation from extra dimensions as an alternative 
explanation for the recent observational data (Deffayet, 
Dvah & Gabadadze, 2002; Alcaniz, 2002; Jain, Dev & 
Alcaniz, 2002; Alcaniz, Jain & Dev, 2002; Lue & Stark- 
man, 2003; Lue, Scoccimarro & Starkman, 2004; Al- 
caniz & Pires, 2004; Zhu & Alcaniz, 2004; Sahni & 
Shtanov 2004). For example, in Sahni & Shtanov (2002) 
a new class of braneworld models which admit a wider 
range of possibilities for dark energy than do the usual 
quintessence scenarios was investigated. For a subclass of 
the parameter values, the acceleration of the universe in 
this class of models can be a transient phenomena which 
could help reconcile an accelerating universe with the re- 
quirements of string/M-theory (Fischler et al. 2001). 

Recently, Dvali & Turner (2003) explored the phe- 
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nomcnology and detectability of a correction on the 
Friedmann equation of the form (1 — ilm)H°' / H"^'^ , 
where fi^ is the matter density parameter, H is the Hub- 
ble parameter (the subscript "o" refers to present time) 
and a is a parameter to be adjusted by the observational 
data. Such a correction behaves like a dark energy with 
an effective equation of state given by a; = — 1 + a/2 in 
the recent past and like a cosmological constant (u) = —1) 
in the distant future. Also based on extra dimensions 
physics, Freese & Lewis (2002) proposed the so-called 
Cardassian expansion, a model in which the universe 
is flat, matter dominated and currently accelerated. In 
the Cardassian universe, the new Friedmann equation 
is given by = g{Pm)i where g{pm) is an arbitrary 
function of the matter energy density pm- In the first 
version of these scenarios, the function g{pm) was given 
by g{pm) = Apm + Bp'^'^^, with the second term driving 
the acceleration of the universe at a late epoch after it 
becomes dominant (a detailed discussion for the origin 
of this Cardassian term from extra dimensions physics 
is given by Chung & Freese 2002. See also Cline & 
Vinet 2003). Although being completely different from 
the physical viewpoint, it was promptly realized that by 
identifying some free parameters Cardassian models and 
quintessence scenarios parameterized by an equation of 
state p = Lup predict the same observational effects in 
what concerns tests involving only the evolution of the 
Hubble parameter with redshift (Freese & Lewis 2002). 

More recently, several forms for the function g{pm) 
have been proposed (Freese 2003). In particular, Wang 
et al. (2003) have studied some observational character- 
istics of a direct generalization of the original Cardas- 
sian model. According to the authors, the observational 
expressions in this new scenario are very different from 
generic quintessence cosmologies and fully determined by 
two dimensionless parameters, n and q. Observational 
constraints from a variety of astronomical data have also 
been investigated recently, both in the original Cardas- 
sian model (Zhu & Fujimoto 2002; 2003; Sen & Sen 
2003a; 2003b) and in its generalized versions (Wang et al. 
2003; Savage, Sugiyama & Freese 2004; Amarzguioui, El- 
garoy & Multamaki, 2004). Perhaps the most interesting 
feature of Cardassian models is that although being mat- 
ter dominated, they may be accelerating and may still 
reconcile the indications for a flat universe (ritotai = 1) 
from CMB observations with clustering estimates that 
point consistently to ri„i — 0.3 with no need to invoke 
either a new dark component or a curvature term. In 
these scenarios, it happens through a redefinition of the 
value of the critical density (Freese & Lewis 2002) - see 
also below. 

The aim of this papcT is to test the viability of gener- 
alized Cardassian (GC) scenarios from statistical prop- 
erties of gravitational lensing. To this end, we use the 
recent Cosmic All Sky Survey (CLASS) lensing sample, 
which constitutes the so far largest lensing sample suit- 
able for statistical analysis. We show that for a large 
interval of the parameters n and q characterizing the 
models, GC scenarios are fully compatible with the cur- 
rent gravitational lensing observations. We also explore 
other observational aspects of the model like, e.g., those 
related to the acceleration redshift (i.e., the redshift at 
which the universe begins to accelerate) and to the age- 
redshift relation. 



2. THE GENERALIZED CARDASSIAN SCENARIO 

The modified Friedmann equation for GC models is 
(Wang et al. 2003) 
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where n and q are two free parameters to be adjusted 
by the observational data, Pcard = Po(l + Zcard)^ is the 
energy density at which the two terms inside the bracket 
become equal and po is the present day matter density. 
Note that for values of n = and g = 1, GC scenarios 
reduce to Cold Dark Matter models with a cosmological 
constant (ACDM). As explained in Wang et al. (2003), 
the first term inside the bracket dominates initially in 
such a way that the standard picture of the early universe 
is completely maintained. At a redshift Zcard, these two 
terms become equal and, afterwards, the second term 
dominates, leading or not to an accelerated expansion 
(note, however, that Zcard is not necessarily the accel- 
eration redshift. See, e.g., the discussion below on the 
deceleration parameter). 
Evaluating Eq. (1) for present day quantities we find 
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From this equation, we see that po is also the critical 
density that now can be written as 
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where Pc.oid = 1-88 x is the standard 

critical density and h is the present day Hubble param- 
eter in units of 100 km.s~^Mpc~^. Note that for some 
combinations of the parameters Zcard, n and q the crit- 
ical density can be much lower than the one previously 
estimated. In other words it means that in the context 
of GC models it is possible to make the dynamical esti- 
mates of the quantity of matter that consistently point to 
Po — (0.2 — 0.4)pc,oid compatible with the observational 
evidence for a flat universe from CMB observations and 
the inflationary flatness prediction with no need of a dark 
energy component (see Freese & Lewis (2002) for a more 
detailed discussion). In Fig. 1 we show a generalized 
version of the Figure 1 of Freese & Lewis (2002) in which 
the plane Zcard — n is displayed for selected values of 
q. The contours are labeled indicating the fraction of 
the standard critical density for different combinations 
of 

Zcard and n. In particular, the points inside the shad- 
owed area delimited by the contours 0.2 - 0.4 are roughly 
consistent with the present clustering estimates (Calberg 
et al. 1996; Dekel et al. 1997). 

From Eq. (3), we see that the observed matter density 
parameter in GC models can be written as 
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For the GC expansion parameterized by n and q, the 
deceleration parameter as a function of the redshift has 
the following form 
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Fig. 1. — Zcard ^ diagrams for the ratio po/Pc old (Eq. 4) and 
selected values of q. The contours arc labeled indicating the corre- 
sponding fraction of the standard critical density. Points inside the 
shadowed area are roughly consistent with the present clustering 
estimates (Calberg et al, 1996). 



Fig. 2. — Deceleration parameter as a function of redshift for 
some selected values of n and q and = 0.3. The horizontal 

line labeled (a) / (d) (qo = 0) divides models with a accelerating 
(a) or decelerating (d) expansion at a given redshift. As discussed 
in the text, the present value of qo is basically determined by the 
value of n. 



where an overdot denotes derivative with respect to time, 
R{t) is the cosmological scale factor and f(z,Vl°^^^,q^n) 
is given by Eq. (7). The behavior of the deceleration pa- 
rameter as a function of redshift for some selected values 
of n and q and ^^'^ = 0.3 is shown in Figure 2. As dis- 
cussed earlier, although completely dominated by mat- 
ter, GC scenarios allow periods of accelerated expansion 
for some combinations of the parameters n and q. Note 
that the present acceleration is basically determined by 
the value of n and that the smaller its value the more 
accelerated is the present expansion for a given value of 
q. For example, for q = 1.5 and n = 0, GC models ac- 
celerate presently faster [qo ~ —0.75) than flat ACDM 
scenarios with JIa = 0.7 {qo — —0.5) although the accel- 
eration redshift is almost identical {za — 0.7) while for 
the same value of q and n = 0.5 we find Qo — —0.13 and 
Za — 0.52. In order to make clear the difference between 
Zcard and Za, for the latter values of q and n, we find di- 
rectly from Eq.(4) Zcard = 1.06. It means that although 
becoming dominant at Zcard — 1 the second term of Eq. 
(1) will drive an accelerated expansion only ^ 1.7 Gyr 
later, at Za — 0.52. 

From the above equations, it also is straightforward to 
show that the agc-rcdshift relation is now given by (the 
total expanding age of the Universe to is obtained by 
taking z = 0) 
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where x' = = {1 + z) ^ is a convenient integration 
variable and the dimensionless function f{x,fl^^,q,n), 



obtained from Eqs. (1) and (4), is written as 
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An interesting study on the total expanding age of the 
Universe in generalized Cardassian scenarios was recently 
presented by Savage ct al. (2004). In light of first year 
WMAP data, it was shown that for a subset of Cardas- 
sian models the age; of the universe varies in the interval 
13.4 - 13.8 Gyr, which is surprisingly close to the stan- 
dard ACDM prediction (13.7 ± 0.2 Gyr). 

The angular diameter distance to a light source, de- 
fined as the ratio of the source diameter to its angular di- 
ameter, is an important concept in lensing statistics. For 
the class of GC models here investigated, the angular di- 
ameter distance, Dls{zl,zs) = Ron{zL,zs)/{l + zs), 
between two objects, for example a lens at Zl and a 
source (galaxy) at zs, reads 
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3. LENSING CONSTRAINTS USING CLASS SAMPLE 

In this Section we present our study of statitical prop- 
erties of gravitational lenses in GC scenarios based on 
the final CLASS well-defined statistical sample (for sim- 
ilar studies in quintessence and quartcssence scenarios 
see, e.g., Chae et al. 2004; Jain et al. 2004). 

The CLASS sample consists of 8958 radio sources out 
of which 13 sources are multiply imaged. Here we work 
only with those multiply imaged sources whose image- 
splittings are known (or likely) to be caused by single 
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galaxies. There are 9 such radio sources: 0218+357, 
0445+123, 0631+519, 0712+472, 0850+054, 1152+199, 
1422+231, 1933+503, 2319+051. We, therefore, follow- 
ing Chae et al. (2002) and Dev, Jain & Mahajan (2004), 
work with a total of 8954 radio source. The sources 
probed by CLASS a.t v = 5 GHz are well represented 
by power-law differential number-flux density relation: 
\dN/dS\ oc (S/S^)"^ with 1] = 2.07 + 0.02 (1.97 + 0.14) 
for 5 > 50 (< 5°) where 5° = 30 mJy (Browne et 
al., 2003). The CLASS unlensed sources can be ade- 
quately described by a Gaussian model with mean red- 
shift z = 1.27 and a dispersion of 0.95. 

In our analysis we assume the singular isothermal 
sphere (SIS) model for the lens mass distribution 
(Turner, Ostriker & Gott 1984). For the present anal- 
ysis we also ignore the evolution of the number density 
of galaxies and assume that the comoving number den- 
sity is conserved and currently given by 



(l>{L)dL, 
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where 4>(L) is the well known Schechter Luminosity Func- 
tion (LF). The differential optical depth of lensing in 
traversing dz^ with angular separation between and 
(j) + d(j) is (Fukugita, Futamase & Kasai, 1990, Turner, 
1990; Fukugita et al., 1992) 
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where the function F* is defined as 
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In Eq. (10) DoL) Dqs and D^s are, respectively, the 
angular diameter distances from the observer to the lens, 
from the observer to the source and between the lens and 
the source (see Eq. 8). In order to relate the characteris- 
tic luminosity to the characteristic velocity dispersion 

, we use the Faber- Jackson relation (Faber & Jackson, 
1976) for early-type galaxies (L, cx w,''), with 7 = 4. For 
the analysis presented here we neglect the contribution 
of spirals as lenses because their velocity dispersion is 
small when compared to ellipticals. 

Two large-scale galaxy surveys, namely, the 2dF 
Galaxy Redshift-Survey (2dFGRS)i and the Sloan Digi- 
tal Sky Survey (SDSS)^ have produced converging results 
on the total LF. These surveys determined the Schechter 
parameters for galaxies (all types) at z < 0.2. For 
our analysis here, we adopt the normalization corrected 
Schechter parameters of the 2dFGRS survey (Folkes et al. 
1999; Chae 2002): a = -0.74,0* = 0.82 x lO-^/jSMpc-^^ 
V* = 185km/s and F* = 0.014. 

The total optical depth t{z) along the line of sight from 
an observer at 2 = to a source at zs is given by 

T{zs) = ^[Dos{l + zs)fRl. (12) 

^ http;//msowww. anu.edu.au/2dFGRS/ 
^ http://www.sdss.org/ 
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Fig. ,3. — Normalized optical depth (t/F*) as a function of the 
source redshift zs for some selected valuesof n and q and a fixed 
value of QS?" = 0.3. 



where Dqs is given by Eq. (8). Figure 3 shows, for the 
fixed value of f2^* = 0.3, the normalized optical depth 
as a function of the source redshift for values of q = 
0.3, 1.0, 2.0 and 3.0 and n = -1.0, -0.5, 0.0 and 0.5. 
Note that a decrease in the value of n at fixed 17^* and 
q tends to increase the optical depth for lensing. For 
q = 1.0 at Zs = 3.0, the value of t/F* for n = 0.5 is 
down from that one for n = —1.0 by a factor of 2.02, 
while the same values of n and q = 2.0 provide values 
for t/F* that are up from the previous ones only by a 
factor of ~ 1.1 and 1.25, respectively. It clearly shows 
that the optical depth is a more sensitive function to the 
parameter n than to the index q. As commented earlier, 
this particular feature is also noted for the other analyses 
discussed in this paper. 
The normalized image angular separation distribution 

,2 

for a source at zs is obtained by integrating g^j^ over 
ZL, i.e., 

— = r ^^dzL. (13) 

d(j) t{zs) Jo dzLd4> 



The corrected (for magnification and selection effects) 
image separation distribution function for a single source 
at redshift zs is given by (Kochanek, 1996; Chiba & 
Yoshii, 2001) 
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Similarly, the corrected lensing probability for a given 
source at redshift z can be written as 

P' = t{zs) J deb, (15) 

where (h and AO are related to as = „ and B 

is the magnification bias. It is worth emphasizing that 
this is considered because, as widely known, gravitational 
lensing causes a magnification of images and this trans- 
fers the lensed sources to higher flux density bins. In 
other words, the lensed sources are over- represented in a 
flux-limited sample. 

The magnification bias B{zs, S^,) increases the lensing 
probability significantly in a bin of total flux density (S^) 
by a factor 

-1 
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Here iV2<,(> 3,^) is the intrinsic flux density relation for 
the source population at redshift zs- N^^ {> S„) gives the 
number of sources at redshift zs having flux greater than 
5*1,. For the SIS model, the magnification probability dis- 
tribution isp(/i) — 8//i'^. The minimum (/imm) and max- 
imum (fJ-max) total magnifications in Eq. (16) depend on 
the observational characteristics as well as on the lens 
model. For the SIS model, the minimum total magnifi- 
cation is fJ-min — 2 while fimax = oo. The magnification 
bias B depends on the differential number-flux density 
relation \dNzc:{> S^)/dS^\, which means that such a re- 
lation needs to be known as a function of the source red- 
shift. Since, at present, redshifts of only a few CLASS 
sources are known, we ignore redshift dependence of the 
differential number-flux density relation. Following Chac 
(2002), we further ignore the dependence of the differen- 
tial number-flux density relation on the spectral index of 
the source. 

Finally, we emphasize that two important selection cri- 
teria for CLASS statistical sample are (i) that the ratio 
of the flux densities of the fainter to the brighter images 



is > 0.1. Given such an observational limit, the 



minimum total magnification for double imaging for the 
adopted model of the lens is /i,„m = 2(1 -I- TZmin/^ — 
T^min)', (ii) that the image components in lens systems 
must have separations > 0.3 arcsec. We incorporate 
this selection criterion by setting the lower limit of A6 
in equation ifT^ as 0.3 arcsec. 

We perform a maximum-likelihood analysis to deter- 
mine the confidence levels in the n — q space. The likeli- 
hood function is defined by 



C^Hil-Pl) l[P'{A0), 



(17) 



fe=i 



where Nl is the observed number of multiple-imaged 
lensed radio sources and Nu is the number of unlensed 
sources in the adopted sample. Pj is the probability of 
the i*'^ source to get lensed and Pj,[A9) is the probability 
of the k*-^ source to get lensed with the observed image 
separation AO. 

Figure 4 shows the main results of our analysis. There, 
we show contours of constant likelihood (68.3% and 
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Fig. 4. — Confidence regions in the plane n — q arising from 
lensing statistics for a fixed value of = 0.3. Solid (dashed) 

lines indicate contours of constant likelihood at 68.3% c.l. (95.4% 
C.I.). The best-fit model is indicated by "x". 



95.4%) in the parameter space n — q. From the above 
equation we find that the maximum value of the like- 
lihood function is located at n = —2.32 and q ~ 0.05, 
which means that the cardassian term of Eq. (1) takes 
over only recently, i.e., at Zcard — 0(10"^). At the la 
level, however, the entire range of q (if we consider for 
instance < g < 2) is compatible with the observa- 
tional data for a fixed value of ri°j^ = 0.3 (a similar 
conclusion also holds for small variations on the value of 
^m'* = 0.3, i.e., ±0.1, which is inside the uncertainties of 
current clustering estimates). As observed earlier, this 
result suggests that a large class of GCs cenarios is in 
accordance with the current gravitational lensing data. 
This best-fit scenario also corresponds to a 7.7/i~^-Gyr- 
old universe {z — 0) and provides an universe old enough 
to accommodate some recent age estimates of high-z ob- 
jects. For example, at z = 1.55 and z — 1.43 Eq. (6) 
provides^ tz — 3.5 Gyr and tz — 3.7 Gyr, respectively, 
i.e., values that are in agreement with the age estimates 
for the radio galaxies LBDS 53W091 and LBDS 53W069 
(Dunlop et al. 1996; Dunlop 1988; Alcaniz & Lima 1999). 
For the recent discovery of the quasar APM 08279-f 5255 
(Hasinger & Komossa 2002; Komossa & Hasinger 2003) 
at z = 3.91, however, these values of q and n provide 
tz — 1.4 Gyr while the age estimate for this object lies 
between 2.0 - 3.0 Gyr. A similar problem is also faced 
by several cosmological scenarios, including the current 
concordance ACDM model (see Friaga, Alcaniz & Lima, 
2004). 

^ For Ho = 64 kms^'^Mpc"^ , which corresponds to Icr lower 
bound obtained by the HST key project (Freedman et al. 2003). 
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4. CONCLUSION 

The possibility of an accelerating universe from dis- 
tance measurements of type la supernovae constitutes 
one of the most important results of modern cosmology. 
These observations naturally lead to the idea of a domi- 
nant dark energy component with negative pressure since 
all known types of matter with positive pressure gener- 
ate attractive forces and decelerate the expansion of the 
universe [R oc —{p + 3p)]. On the other hand, the re- 
alization that dark energy or the effects of dark energy 
could be a manifestation of a modification to the Fried- 
mann equation arising from extra dimension physics has 
opened up an unprecedented opportunity to establish a 
more solid connection between particle physics and cos- 
mology. Many braneworld scenarios have been proposed 
in the recent literature with some of them presenting in- 
teresting features which make them a natural alternative 
to the standard model. Here we have analyzed some ob- 
servational consequences of one of these scenarios, the 
so-called generalized cardassian expansion, recently pro- 
posed by Wang et al. (2003). We have studied the ob- 
servational constraints on the parameters n and q that 
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fully characterize the model from statistical properties 
of gravitational lensing. To this end, we have used the 
final CLASS well-defined statistical sample, which con- 
stitutes the most recent gravitational lensing data. From 
our analysis we have found that at la level a large por- 
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